function f_quadratic = cal_fq_Closed_Form(N,K,uu,vv,ww,hGG,noise)
f_quadratic=0;
for nn=1:N
    for kk=1:K
        sum_wmlH_hmnk = 0;
        for mm=1:N
            for ll=1:K
                sum_wmlH_hmnk = sum_wmlH_hmnk + abs(ww(:,mm,ll)'*hGG{mm,nn,kk})^2;
            end
        end
        f_quadratic=f_quadratic + 2 * sqrt((1+uu(nn,kk))) * real(ww(:,nn,kk)'*hGG{nn,nn,kk}*vv(nn,kk)) ...
            - (vv(nn,kk)' * vv(nn,kk)) * (noise + sum_wmlH_hmnk);
    end
end
end